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Abstract 

The B^-B^ oscillation frequency Am^ is measured with 36 pb~^ of data collected in 
pp collisions at ^/s = 7 TeV by the LHCb experiment at the Large Hadron Collider. 
A total of 1381 Bg^ D~tt^ and -B^ — >• D~ tt~^ Tr~ tt'^ signal decays are reconstructed, 
with average decay time resolutions of 44 fs and 36 fs, respectively. An oscillation 
signal with a statistical significance of 4.6 o" is observed. The measured oscillation 
frequency is Arus = 17.63 ± 0.11 (stat) it 0.02 (syst) ps~^. 
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1 Introduction 



After the observation of B^-B^ mixing and the measurement of its strength in 1987 [1], it 
took a further 19 years for the B^-B^^ frequency to be measured for the first time [2], [3]. 
This is mainly due to the fact that the B^-B^ oscillation frequency is 35 times larger than 
that for the B^-B^ system, posing a considerable challenge for the decay time resolution 
of detectors. For the LHCb experiment, the abihty to resolve these fast -B°-S° oscillations 
is a prerequisite for many physics analyses. In particular it is essential for the study of 
the time-dependent CP asymmetry of B^ — )■ J/ip decays [4] . The oscillation frequency 
in the B^-B^ system is given by the mass difference between the heavy and light mass 
eigenstates, Am^ (we use units with ^=1). In this letter, we report a measurement of 
ArUs by the LHCb experiment with data collected in 2010. 

The LHCb spectrometer covers the pseudo-rapidity range 2 to 5. In this region, b 
hadrons are produced with a large Lorentz boost and have an average flight path of 
7 mm. The LHCb detector consists of several components arranged along the LHC 
beam line. The vertex detector (VELO) surrounds the collision point, followed by a flrst 
Ring Imaging Cherenkov (RICH) counter, a tracking station, a dipole magnet, three more 
tracking stations, a second RICH detector, a calorimeter system and a muon detector. The 
calorimeter system consists of a scintillating pad detector (SPD), a preshower detector, 
an electromagnetic calorimeter and a hadronic calorimeter. A detailed description of the 
detector can be found in Ref. [5]. The precise spatial resolution of the VELO results 
in an impact parameter resolution of 20-50 /xm in the x and y directions^ for charged 
particles with transverse momenta in the range relevant for B^ daughter tracks used in 
this analysis. The x and y resolution in the position of the primary vertex reconstruction 
is about 15 fim while the z resolution is about 80 /xm. This excellent performance results 
in the decay time resolution needed to observe the fast B^-B'^ oscillations. The invariant 
mass resolution provided by the tracking system and the tt/K separation given by the 
two RICH detectors provide clean 5° meson signals with small background. The particle 
identiflcation capabilities of the RICH together with the calorimeter and muon systems 
allow the initial flavour of the B^ to be tagged using charged kaons, electrons and muons, 
respectively. 

In the next section, the data sample used and the analysis strategy are introduced. 
This is followed by descriptions of the analysis of the invariant mass and decay time 
distributions, and the flavour tagging. Finally, we discuss the fit result for the oscillation 
frequency and the associated systematic uncertainties. 

"'^LHCb uses a right-handed Cartesian coordinate system with the x direction pointing inside the LHC 
ring, the y direction pointing upwards and the z direction running along the beamhne from the interaction 
point towards the spectrometer. 
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2 Data sample and analysis strategy 



The analysis uses candidates reconstructed in four flavour-specific decay modes, ^ 
namely 5° ^ Dj{(t){K+K-)7r-)7i+, 5° ^ Dj {K*%K+7r-)K-)7r+ , 
D~{K~^K~Tr~)Tr~^ and B^ D~ {K^ K~t:~)t:^t:~t:^ . To avoid double counting, can- 
didates that pass the selection criteria of one mode are not considered for the following 
modes. All reconstructed decays are flavour-specific final states, thus the flavour of the 
B^ at the time of its decay is given by the charges of the final state particles of the decay. 
A combination of tagging algorithms is used to identify the 5° flavour at production. 
The algorithms provide for each event a tagging decision as well as an estimate of the 
probability that this decision is wrong (mistag probability). These algorithms have been 
optimized and calibrated using large event samples of flavour-specific B — > jji^D*~X and 
S"^— >■ J/il^K^ decays and a sample of > D^n^ decays. 

The analysis is based on a data set of 36 pb^^ of pp collisions at ^/s = 7 TeV collected 
in 2010. The first trigger level is implemented in hardware, while the second trigger level is 
based on software. Trigger conditions were progressively tightened over the duration of the 
data taking period to cope with the rapidly increasing instantaneous luminosities delivered 
by the LHC. In the hardware trigger, the events used in this analysis were selected by 
requiring a cluster with a minimum transverse energy in the hadronic calorimeter. The 
applied threshold was increased from 2.5 to 3.6 GeV throughout the data taking period. 
A cut on the number of hits in the SPD detector was applied to reject very high occupancy 
events. The software trigger for the first 2.4 pb~^ of data required a good quality displaced 
vertex reconstructed from two tracks with transverse momenta pt of at least 500 MeV/ c. 
For the remaining data, a two-level software trigger was applied. A good quality track 
with large impact parameter with respect to the primary vertex was required with p^ > 
1.85 GeV/c and momentum p > 13.3 GeV/c [6]. For events passing these criteria, a good 
quality displaced vertex was required, formed out of two tracks with px > 0.5 GeV/c and 
p > 5 GeV/c and with a mass variable in the range 2 to 7 GeV/c^ [7]. 

Some of the offline event selection criteria are optimized individually for each of the 
four decay modes under study. In this way specific features such as the masses of the 
intermediate (p and K*^ resonances or the helicity angle distribution of the K*^ can be 
used. The selection criteria common to all decay modes exploit the long B^ hfetime 
by applying cuts on the impact parameters of the daughter tracks, on the angle of the 
reconstructed B^ momentum relative to the line between the reconstructed primary vertex 
and the B^ vertex and on the B^ decay time. Additional cuts are applied on the p and pt 
of the Bg candidate and its decay products as well as on particle identification variables 
and on track and vertex quality. Finally, cuts on the impact parameter significance of the 
reconstructed Dj and its distance of closest approach to the primary vertex are applied. 
The reconstructed D~ mass is required to be consistent with the PDG value [8]. After 
this selection, a total of about 14,400 candidates remain in the B^ — )> Djtt^ invariant 
mass window of [4.80, 5.85] GeV/c^ and in the B^^ D~TT~^7i^n^ invariant mass window 
of [5.00, 5.60] GeV/c^. 

^Unless explicitly stated, inclusion of charge-conjugated modes is implied. 
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An unbinned likelihood method is employed to fit simultaneously the invariant mass 
and decay time distributions of the four decay modes. The probability density functions 
(PDFs) for the signal and for the background in each of the four modes can be written as 

V = Vm{m) Vtit, q\at, rj) Vr^iv), (1) 

where m is the reconstructed invariant mass of the candidate, t is its reconstructed 
decay time and at is the event-by-event estimate of the decay time resolution given by the 
event reconstruction algorithm. The tagging decision q can be (no tag), —1 (different 
flavour at production and decay) or -|-1 (same flavour at production and decay). The 
predicted event-by-event mistag probability r) can take values between and 0.5. The 
terms Vm and Vt describe the invariant mass distribution and the decay time distribution, 
respectively. Vt is a conditional probability depending on cr^ and f]. The terms V^t and 
are required to ensure the proper relative normalization of Vt for signal and background [9] . 
These terms are determined directly from the data, using the measured distribution in the 
upper Bg invariant mass sideband for the background PDF and the sideband subtracted 
distribution in the invariant mass signal region for the signal PDF. 

3 Fit to the invariant mass distributions 

The invariant mass of each candidate is determined in a vertex fit using a constraint 
on the D~ mass. The invariant mass spectra for the four decay modes after all selection 
criteria are shown in Fig. 1. The four distributions are fit simultaneously taking into 
account contributions from signal, combinatorial background and b decay backgrounds. 
The signals are described by Gaussian distributions. The fit constrains the mean of the 
Gaussian distributions to be the same for all four decay modes, whereas it allows the width 
to be different for the B^ Djn^ and the B^ — )■ Djn'^n'n'^ modes, respectively. The 
combinatorial backgrounds are described by exponential functions. Their parameters are 
allowed to vary individually for the four decay modes. An alternative parameterization of 
the combinatorial backgrounds by a first order polynomial is used as part of the systematic 
studies. 

The b decay backgrounds include partially reconstructed B^ decays, as well as fully 
and partially reconstructed S° and A;, decays with one mis-identified daughter particle. 
Their shapes are derived from a large simulated event sample, where all selection cuts were 
applied on generator level quantities. The invariant mass spectra were then smeared with 
a Gaussian distribution to take into account effects of detector resolution. This approach 
was validated by comparing the results with those from a full simulation including a 
detailed description of the detector response. The relative normahzation factors for the 
different b decay backgrounds are parameters in the fit. They are constrained to be the 
same for the three -B^— ?■ D~tt^ decay modes. 

The fit returns a value of m(5°) = 5364.7 ± 0.7 MeV/c^, about 1.5 McV/c^ below the 
PDG value [8]. This mass shift is attributed to imperfections in the detector alignment 
and magnetic field calibration. A dedicated study on the momentum scale resulted in 
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Figure 1: Mass distributions for a) S° D;{(^n-)n+, b) {K*'^K-)7i+ , 

c) Bg-^ Dj{K~^K~tt~)tt'^ and d) Dj-n-^n^'n-^ candidates. The fits and the various 

background components are described in the text. "Partial" refers to background from 
partially reconstructed B^ decays, "mis-id" refers to background from fully or partially 
reconstructed B^ and A;, decays with one mis-identified daughter particle, and "comb" 
refers to combinatorial background. 



a correction for this effect [10]. This calibration procedure is however not used for the 
analysis presented here as the momentum scale correction largely cancels in the calcula- 
tion of ArUs. The mass templates describing b decay backgrounds are shifted according 
to the observed bias. The fit gives signal mass resolutions of am = 18.1 MeV/c^ for the 
5°— 7- Djn^ modes and o"m = 12.7 MeV/c^ for the >■ D~7i~^n^7r^ mode, respectively. 
The signal yields extracted from the fit are summarized in Table 1. For the remainder 
of the analysis, the invariant -B^ mass range is limited to [m{Bg) — 3am, 5.85 GeV/c^] 
and [m{B{) - 3ct^, 5.60 GeV/c^] for the 5° ^ D-tt+ and 5° -> D;7r+TT-n+ modes, 
respectively. The lower cut of this asymmetric mass window is chosen to reject all back- 
ground candidates from partial reconstructed B^ decays. The only remaining b decay 
backgrounds are thus due to mis-identified B^ and Ab decays. The candidates in the high 
mass sidebands provide a clean sample of combinatorial background. Including them in 
the fit permits to determine the decay time distribution and tagging behaviour of this 
background contribution. 
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Table 1: S° signal yields. 



Decay mode 


Signal yield 


5°^ D;(07r-)7r+ 
5°^ D-{k*^K-)ti+ 
D-{K+K-7r-)n+ 


515 ± 25 
338 ± 27 
283 ± 27 
245 ± 46 


Total 


1381 ± 65 



The parameters derived in the fit to the mass distributions are fixed for the remainder of 
the analysis. 

4 Fit to the decay time distribution 

Ignoring detector resolution effects, selection biases and flavour tagging, the distribution 
of the decay time t of the signal is described by 

Vtit) oc r,e-r»* cosh (2) 

where Tg is the decay width and AF^ the decay width difference between the heavy 
and the hght mass eigenstates. In the fit AF^ is fixed to its PDG value of 0.09 F^ [8]. As 

part of the evaluation of systematic uncertainties on Arris, the assumed vahie of AF^ is 
varied within its current uncertainty between and 0.2 F^. The step function 6{t) restricts 
the PDF to positive decay times. 

The true decay time is convolved with the decay time resolution function of the detec- 
tor. An event-by-event estimate of the decay time resolution is calculated by the fitting 
algorithm, which reconstructs the decay vertex of the B^ and computes its decay length 
and decay time. No constraint on the D~ mass is applied in the computation of the 
decay time in order to minimize sensitivity to the knowledge of the momentum scale of 
the experiment. The decay time uncertainty calculated by the fitting algorithm does not 
include possible effects from an imperfect understanding of the detector material or its 
spatial alignment. To correct for such effects, the calculated event-by-event decay time 
uncertainties, at, are multiplied by a constant scale factor S^f The value of S^t is de- 
termined from data, using a sample of fake B^ candidates formed by a prompt Dj and 
a 7r+ from the primary vertex. The contamination due to secondary Dj from B decays 
is estimated and statistically subtracted using the measured Dj impact parameter dis- 
tribution. The distribution of decay times for this fake 5^ sample, each divided by its 
calculated event-by-event uncertainty, is fitted with a Gaussian function and is taken 
as the resulting standard deviation. Using the full sample of fake B^ candidates, a value 
of 5*^4 = 1.3 is obtained. This value is used as the nominal scale factor in the Anig analysis. 
Studying different regions of phase space of the fake B^ candidates separately, values for 
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Sfjt between 1.2 and 1.4 are obtained. This variation is taken into account for evaluat- 
ing the systematic uncertainties on Am^. Including the nominal scale factor 
the average decay time resolution is 44 fs for the — > D'n'^ sample and 36 fs for the 
Bg — >■ Dj7r'^n~7r^ sample. The decay time resolution is taken into account in the PDF 
by convolving Eq. 2 with a Gaussian G with mean zero and standard deviation 1.3 at. 

The shape of the decay time distribution is distorted by trigger and offline selection 
criteria which require several particles with large impact parameter with respect to the 
primary vertex. This is accounted for in the PDF by introducing an acceptance function 
e{t), derived from a full detector simulation. Determining e(t) from simulation is deemed 
acceptable since it cancels to first order in the determination of Am^. The untagged signal 
decay time PDF becomes 

^G{t,Sa,at)e{t). (3) 

The decay time distributions for the b decay backgrounds from and decays are 
described in the same way as that for signal B^ candidates, using the PDG values for 
their lifetimes and Ar=0. The shape of the decay time distribution for the combinatorial 
background is described by the sum of two exponential functions multiplied by a second 
order polynomial. The parameters of these functions are derived from the high mass side- 
bands. Figure 2 illustrates the results of the lifetime fit. Within its statistical uncertainty 
the reconstructed B^ lifetime agrees with the PDG value [8]. 

5 Flavour tagging 

To determine the fiavour of the B^ candidate at production we exploit the fact that b 
quarks are predominantly produced in quark-antiquark pairs. The quark which is not 
part of the B^ meson gives rise to an opposite-side b hadron. For opposite-side b hadron 
decay candidates, the charge of displaced muons, electrons and kaons and a decay vertex 
charge estimate arc combined using a neural network to form a single opposite-side tagging 
decision. The tagging decision has a probability to be wrong which is called the mistag 
probability, ou. For each event an estimate, rj, of the mistag probability, is determined 
based upon topological and kinematic properties of the event, including the number of 
primary vertices, the number of tagging particle candidates, the impact parameter of the 
tagging particle and of the B^ candidate with respect to the primary vertex, and the p and 
Pt of the selected tagging particle and the B^ candidate. The optimization of the tagging 
algorithms and an initial calibration of rj are performed in an independent analysis using 
large event samples of S — > ^^D*~X and B^ — >■ J/ipK'^ decays. More details on the 
individual tagging algorithms and this calibration procedure can be found in Ref. [11]. 

The B — 7> fi^D*~X and B^ J/ipK^ events used in the optimization and calibration 
were collected using different trigger and selection criteria than for the B^ Djn^ and 
Bg L>~7r+7r~7r+ events used in the Am^ analysis described here. As trigger and selection 
cuts can bias the distributions of the event properties used by the tagging algorithms. 



Vt{t\at) (X 
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decay time [ps] decay time [ps] 

Figure 2: Decay time distributions for a) 5°^ {(j)Tx-)Ti^ , b) {K*^ K-)ti+ , 

c) 5° D-{K+K-tt-)tt+ and d) 5° -> D77r+7r-7r+ candidates. The data and the fit 
projection are from a mass range of ± 3am around the reconstructed mass. The 
abbreviations for the various fit components are introduced in Fig. 1. 



this could resuh in a biased estimate for the — > Djtv^ and B^ — ?► D J 71^11-^71^ events. 
Therefore, a re-cahbration is performed using a sample of 6,000 B^ — )■ D^tt^ events, 
which have a similar topology to the B^ Djtt^ and B^ — )■ Djti^ti^ii^ events, and 
were collected using the same trigger and similar selection cuts. This event sample is used 
to perform a measurement of the B^-B^ flavour oscillation using a very similar method 
to that described here. In that measurement the true event mistag probability, w, is 
parameterized as a linear function of t] using the relationship uj{rj) = a + b x [1] — {r])) , 
where (rj) = 0.3276 is the mean of the distribution of the 77 values obtained from the 
initial tagger optimization. The parameters a = 0.311 ± 0.022 and b = 0.61 ± 0.25 
are determined as part of the maximum likelihood fit of the B^-B^ oscillation signal 
and found to be consistent with the original calibration. As a by-product of this re- 
calibration procedure the B^-B^ oscillation frequency is measured. The resulting value 
of Anid = 0.499 ± 0.032 (stat) ± 0.003 (syst) ps~^, though statistically less precise, is 
in good agreement with the PDG value of Arrid = 0.507 ± 0.004 ps^^ [8] and provides a 
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valuable cross check of the procedure. 

The statistical power of the tagging is determined by the "effective" tagging efficiency 
for signal events and is defined as 

eeff = X - 2u{r^i)f x W,, (4) 

where the signal tagging efficiency eg is a free parameter in the fit of the oscillation 
frequency described in the next section. Wj is the probability for being a signal event as 
determined by the invariant mass and decay time PDFs. The index i runs over all 
candidates. 

6 Measurement of the oscillation frequency 

To determine the oscillation frequency, Amg, the decay time PDF for signal candidates 
with tagging information is modified in the following way: 



cosh ^^^^^ + ? [1 ~ 2c(;(?7)] cos(Amst) 



9{t) 



2 

^ G{t, Sa, at) €{t) (5) 

The decay time PDF for untagged signal events is given by Eq. (3) multiphed by an 
additional factor (1 — eg). The calibration parameters a and b of the mistag probability 
uj{ri) are identical for all signal and b decay background components. Within Gaussian 
constraints they are set to the values found in the calibration described in the previous 
section. The signal tagging efficiency eg for the -B°— >■ -Dj7r+ and -B°— )> D~7r'^7r~7r'^ modes 
are two separate parameters in the fit. The same values of eg are however used for signal 
and b decay background components in each of these two categories. In the description of 
the combinatorial background a separate parameter for the tagging efficiency is introduced 
for each of the four modes. In addition, tagging asymmetry parameters are introduced 
in the PDFs for the combinatorial background, to allow for a different number of events 
tagged as or in each mode. As expected the fit results for these asymmetries are 
compatible with zero. 

The fit for the oscillation frequency Am^ is performed simultaneously to all four 
decay modes and gives Am^ = 17.63 ± 0.11 ps~^ (statistical uncertainty only). Signal 
tagging efficiencies of eg = (23.6 ± 1.3) % and eg = (17.6 ± 3.2) % are found for the 
Bg — > Djn'^ and B^ — > £)j7r+7r~7r+ modes, respectively. The combined effective tagging 
efficiency for all four modes is eeff = (3.8 ± 2.1) %. The likelihood profile as a function of 
the assumed oscillation frequency Arris is shown in Fig. 3. The statistical significance of 
the signal is evaluated to be 4.6 ex by comparing the likelihood value at the minimum of 
the fit with that found in the limit Am^ = oo. 

To illustrate the oscillation pattern, we define the time dependent mixing asymmetry 

as 
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™^^ ^~iV+(t) + iV-(t) ^ ^ 

where N^{t) and N~{t) are the number of background subtracted signal candidates 
with a given decay time t and tagging decision +1 and —1, respectively. Note, that this 
definition of the asymmetry does not include any information on the mistag probabilities 
and therefore does not use the full information of the likelihood fit. Despite the limited 
size of the sample, the oscillation pattern is clearly visible when the asymmetry is plotted 
in bins of the decay time modulo 2tt/ Arris (Fig. 4). In an ideal scenario of perfect tagging 
and perfect decay time resolution the amplitude of this oscillation would be 1.0. The 
observed amplitude is reduced due to the performance of the tagging algorithm by a 
factor 0.41. Another reduction of 0.65 occurs due to the limited decay time resolution. 



7 Systematic uncertainties 

The dominant source of systematic uncertainty is due to the knowledge of the absolute 
decay time scale of the experiment. This uncertainty is dominated by the knowledge of 
the z scale. A relative uncertainty of 0.1% on the z scale and thus on the decay length 
is assigned based on comparisons of detector surveys and a software alignment using re- 
constructed tracks. This leads to a systematic uncertainty of 0.018 ps~^ on Arris. A 
second contribution to the decay time scale is due to the momentum scale of the exper- 
iment. From an independent analysis of the mass scale using various known resonances 
an uncertainty of the uncalibrated momentum scale of less than 0.1% is estimated. This 
uncertainty partially cancels as it enters both the reconstructed B'^ mass and the B^ mo- 
mentum. The resulting relative uncertainty on the decay time is 0.02%, which translates 
to an absolute systematic uncertainty of 0.004 ps~^ on Am^. 



9 




0.1 0.2 0.3 

t modulo 2k I Arrig [ ps ] 



Figure 4: Measured asymmetry for 5° candidates in bins of the decay time t modulo 
27r/Am5. The projection of the hkehhood fit is superimposed. 



Table 2: Summary of the systematic uncertainties on Amg- The total systematic uncer- 
tainty is defined as the quadratic sum of the individual components. 



Source 


Uncertainty [ps ^] 


Momentum scale 


0.004 


z scale 


0.018 


Comb, background mass shape 


0.010 


Decay time resolution 


0.006 


Total systematic uncertainty 


0.022 



The next largest systematic uncertainty is related to the description of the combinato- 
rial background in the fit to the mass spectra. It is evaluated by replacing the exponential 
function by a first order polynomial. Based on the shift in the value obtained for Am^, 
a systematic uncertainty of 0.010 ps~^ is assigned. Finally, based on variations of the 
decay time resolution scale factor S^^^ within its estimated uncertainty from 1.2 to 1.4, 
a systematic uncertainty of 0.006 ps~^ is assigned on Am^. These contributions to the 
systematic uncertainty on Am^ are summarized in Table 2. 

Various other possible sources of systematic effects have been studied, such as the decay 
time resolution model, the decay time acceptance, releasing parameters of the invariant 
mass and decay time PDF in the mixing fit, different parameterizations of the invariant 
mass of the b decay backgrounds and variations of the value of AF^. They are found to 
be negligible. 
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8 Conclusion 



A measurement of the B^-B^ oscillation frequency Am^ is performed using B^ — )■ Dj7r~^ 
and — )■ D~7r~^n~n^ decays collected in 36 pb~^ of pp collisions at ^/s — 7 TeV in 2010. 
The result is found to be 

Arus = 17.63 ± 0.11 (stat) ± 0.02 (syst) ps"\ (7) 

This is in good agreement with the previous best measurement of ArUs — 17.77 ± 
0.10 (stat) ± 0.07 (sys) ps~^, reported by the CDF collaboration [3]. As a by prod- 
uct of the analysis wc also determine a value for the B^-B^ oscillation frequency 
Arrid = 0.499 ± 0.032 (stat) ± 0.003 (syst) ps^^. Our results are completely domi- 
nated by statistical uncertainties and thus significant improvements are expected with 
larger data sets. 
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